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Abstract 

Diffractive dissociation of particles can be used to study their light-cone wave func- 
tion. Results from Fermilab experiment E791 for diffractive dissociation of 500 
GeV/c ir~ mesons into di-jets are presented. The results show that the \ qq) light-cone 
asymptotic wave function describes the data well for Q 2 ~ 10 (GeV/c) 2 or more. 
Evidence for color transparency comes from a measurement of the ^-dependence of 
the yield of the diffractive di-jets. It is proposed to carry out similar studies for the 
light-cone wave function of the photon. 



Introduction 

The Pion Light- Cone Wave Function 

The internal momentum distribution of valence quarks in hadrons are fun- 
damental to QCD. [1]. They are generated from the valence light-cone wave 
functions integrated over k t < Q 2 , where k t is the intrinsic transverse momen- 
tum of the valence constituents and Q 2 , the total momentum transfer squared. 
Even though these amplitudes were calculated about 20 years ago, there have 
been no direct measurements until those reported here. 

The pion wave function can be expanded in terms of Fock states: 

* = a-ilqq) + a 2 \qqg) + a 3 \qqgg) + .... (1) 



Two functions have been proposed to describe the momentum distribution am- 
plitude for the quark and antiquark in the \qq) configuration. The asymptotic 
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function was calculated using perturbative QCD (pQCD) methods [2-4], and 
is the solution to the pQCD evolution equation for very large Q 2 (Q 2 — > oo): 

<j) as (u) = V3u(l-u). (2) 

u is the fraction of the longitudinal momentum of the pion carried by the quark 
in the infinite momentum frame. It should not be confused with xbj which is 
not specific to a single Fock state. The antiquark carries a fraction (1 — u). 
Using QCD sum rules, Chernyak and Zhitnitsky [5] proposed a function that 
is expected to be correct for low Q 2 : 

<p cz (u) = 5y/3u(l -u)(l -2m) 2 . (3) 

As can be seen from eqns. 2 and 3 and from Fig. 1, there is a large difference 
between the two functions. Measurements of form factors are insensitive to 
the wave function as these quantities are derived by integrating over the wave 
function and interpretation of the results is model dependent [6]. In this work 
we describe an experimental study that maps the momentum distribution of 
the valence \qq) in the pion. This provides the first direct measurement of the 
pion light-cone wave function (squared). The concept of the measurement is 
the following: a high energy pion dissociates diffractively on a heavy nuclear 
target. This is a coherent process in which the quark and antiquark break 
apart and hadronize into two jets. If in this fragmentation process the quark 
momentum is transferred to the jet, measurement of the jet momentum gives 
the quark (and antiquark) momentum. Thus: u measure d = - - ■ From 
simple kinematics and assuming that the masses of the jets are small compared 
with the mass of the di-jets, the virtuality and mass-squared of the di-jets are 
given by: Q 2 ~ M^j = ^t^j where k t is the transverse momentum of each jet. 
By studying the momentum distribution for various k t bins, one can observe 
changes in the apparent fractions of asymptotic and Chernyak-Zhitnitsky (CZ) 
contributions to the pion wave function. 

The basic assumption that the momentum carried by the dissociating qq is 
transferred to the di-jets was examined by Monte Carlo (MC) simulations of 
the asymptotic and (CZ) wave functions (squared). The MC samples were al- 
lowed to hadronize through the LUND PYTHIA-JETSET model [7] and then 
passed through simulation of the experimental apparatus (described in the 
next section) to simulate the effect of unmeasured neutrals and other exper- 
imental distortions. In Fig. 1 the initial distributions at the quark level are 
compared with the final distributions of the detected di-jets. As can be seen, 
the qualitative features of the distributions are retained. The results of this 
analysis come from comparing the observed -u-distribution to a combination 
of the distributions shown, as examples, on the right of Figure 1. 
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Fig. 1. Monte Carlo simulations of squares of the two wave functions at the quark 
level (left) and of the reconstructed distributions of di-jets as detected (right). <j>\ sy 
is the asymptotic function (squared) and c/)q Z is the Chernyak-Zhitnitsky function 
(squared). The di-jet mass used in the simulation is 6 GeV/c 2 and the plots are for 
1.5 GeV/c < k t < 2.5 GeV/c. 

The Color Transparency Effect 

The Color- Transparency (CT) phenomenon is derived from the prediction 
that color fields cancel for physically small color singlet systems of quarks 
and gluons [8]. This effect of color neutrality (or color screening) is expected 
to lead to the suppression of initial and final state interactions for a small 
sized system or point-like configuration (PLC) formed in a large angle hard 
process [9]. Observation of CT requires that a PLC is formed and that the 
energies are high enough so that expansion of the PLC does not occur [10] (the 
frozen approximation). Under conditions of k t >1.5 GeV/c, which translates 
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to Q 2 ~ 10 (GeV/c) 2 and (r) ~ 0.1 fm, observation of these effects can be 
expected. Bertsch et al. [4] proposed that the small \qq) component will be 
filtered by the nucleus. Frankfurt et al. [11] show that for k t > l.5GeV/c the 
interaction with the nucleus is completely coherent and a(\qq)N — ^di-jets N) 
is small. This leads to an A 2 dependence of the forward amplitude squared. 



Experimental Results 



Fermilab experiment E791 recorded 2 x 10 10 events from interactions of a 500 
GeV/c 7r~ beam with carbon and platinum targets. Details of the experiment 
are given in [12]. Only about 10% of the E791 data was used for the analysis 
presented here. The data were analysed by selecting events in which 90% of 
the beam momentum was carried by charged particles. Jets were identifies 
using the JADE jet-finding algorithm [13]. To insure clean selection of two-jet 
events, a minimum k t of 1.25 GeV/c was required and their relative azimuthal 
angle, which for pure di-jets should be 180° was required to be within 20° of 
this value. 

Diffractive di-jets were identified through the e~ bq t dependence of their yield 
(q 2 is the square of the transverse momentum transferred to the nucleus and 
b = <R 3 > where R is the nuclear radius). Figure 2 shows the q 2 distributions 
of di-jet events from platinum and carbon. The different slopes in the low q 2 t 
coherent region reflect the different nuclear radii. Events in this region come 
from diffractive dissociation of the pion. 
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Fig. 2. q 2 distributions of di-jets with 1.5 < k t < 2.0 GeV/c for the platinum 
and carbon targets. The lines are fits of the MC simulations to the data: coher- 
ent dissociation (dotted line), incoherent dissociation (dashed line), background 
(dashed-dotted line), and total fit (solid line). 
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The Pion Wave Function 



For measurement of the wave function we used data from the platinum target 
as it has a sharp diffractive distribution and low background. We used events 
with q 2 < 0.015 GeV/c 2 . For these events, the value of u was computed from 
the measured longitudinal momentum of each jet. A background, estimated 
from the u distribution for events with larger q 2 was subtracted. The analysis 
was carried out in two windows of k t : 1.25 GeV/c < k t < 1.5 GeV/c and 
1.5GeV/c < k t < 2.5 GeV/ c. The resulting u distributions are shown in Fig. 
3. In order to get a measure of the correspondence between the experimental 
results and the calculated light-cone wave functions, we fit the results with 
a linear combination of squares of the two wave functions (right side of Fig. 
1). This assumes an incoherent combination of the two wave functions and 
that the evolution of the CZ function is slow (as stated in [5]). The results 




Fig. 3. The u distribution of diffractive di-jets from the platinum target for 
1.25 < k t < 1.5 GeV/c (left) and for 1.5 < k t < 2.5 GeV/c (right). The solid 
line is a fit to a combination of the asymptotic and CZ wave functions. The dashed 
line shows the contribution from the asymptotic function and the dotted line that 
of the CZ function. 

of the fits are given in Fig. 3 and in Table 1 in terms of the coefficients a as 
and a cz representing the contributions of the asymptotic and CZ functions, 
respectively. The results for the higher k t window show that the asymptotic 
wave function describes the data very well. Hence, for k t >1.5 GeV/c, which 
translates to Q 2 ~ 10 (GeV/c) 2 , the pQCD approach that led to construction 
of the asymptotic wave function is reasonable. The distribution in the lower 
window is consistent with a significant contribution from the CZ wave function 
or may indicate contributions due to other non-perturbative effects. 
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1.25 - 1.5 
1.5 - 2.5 


0.64 
1.00 


±0.12 
±0.10 


±0.07 -0.01 
±0.00 -0.10 


±0.14 -0.12 
±0.10 -0.14 


0.36 
0.00 


T0.12 

to.io 


-0.07 ±0.01 
-0.00 ±0.10 


-0.14 ±0.12 
-0.10 ±0.14 



Table 1 



Asymptotic (a as ) and CZ (a cz ) wave functions contributions to a fit of the data. 
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kt GeV/c k, GeV/c 

Fig. 4. Comparison of the kt distribution of acceptance-corrected data with fits to 
cross section dependence (a) according to a power law, (b) based on a nonper- 
turbative Gaussian wave function for low kt and a power law, as expected from 
perturbative calculations, for high kt- 

The k t dependence of diffractive di-jets is another observable that can show 
how well the perturbative calculations describe the data. As shown in [11] 
assuming interaction via two gluon exchange and 4> as would lead to ~ k^ 6 . 
The results, corrected for experimental acceptance, are shown in figure 4(a) 
fitted by k r t l for k t > 1.25 GeV with n = -9.2±0 A(stat) ±0.3(sys) and x 2 /dof 
= 1.0. This slope is significantly larger than expected. However, the region 
above k t ~ 1.8 GeV/c can be fitted (Fig. 4(b)) with n = — 6.5±2.0 with x 2 /dof 
= 0.8, consistent with the predictions. This would support the evaluation of 
the light-cone wave function at large k t as due to one gluon exchange. The 
lower k t region can be fitted with the non-perturbative Gaussian function: 
iP ~ e -^l [14], w ith (3 = 1.78 ± 0.05(stat) ± O.l(sys) and x 2 /dof = 1.1. 
Model-dependent values in the range of 0.9 - 4.0 were used [14]. These results 
are consistent with the measurements of the wave function that indicated 
noticeable non-perturbative effects up to k t ~ 1.5 GeV/c. 
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kt bin 
Lre V /c 


a 




A r\ 

sys 




a (0 1 J 


1.25 - 1.5 


1.64 


±0.05 


±0.04 -0.11 


±0.06 -0.12 


1.25 


1.5 - 2.0 


1.52 


±0.09 


±0.08 


±0.12 


1.45 


2.0 - 2.5 


1.55 


±0.11 


±0.12 


±0.16 


1.60 



Table 2 

The exponent in a oc ^4 Q , experimental results for coherent dissociation and the 
Color- Transparency (CT) predictions. 

The Color Transparency Effect 

To study the CT effect we measure the A-dependence of the diffractive di-jet 
yield. The coherence length is estimated using 2p; a t, = 1000 GeV/c and Mj ~ 
5 GeV/c 2 . The result is l c ~ 10 fm, larger than the platinum nuclear radius. 
In order to correct for experimental acceptance, we generate MC simulations of 
diffractive di-jets using the asymptotic wave function and di-jet masses of 4,5, 
and 6 GeV/c. The simulated coherent q\ distributions of the di-jets represent 
the nuclear form factors of carbon (R=2.44 fm) and platinum (R=5.27 fm) [15] 
and the incoherent dissociation is simulated according to the nucleon radius 
(R=0.8 fm) [15] truncated at q% < 0.015. A combination of these simulations 
was used to fit the data (Fig. 2). We derive the numbers of producted di-jet 
events in the data for each target in three kt bins by integrating over the 
diffractive terms in the fits. Using the resulting yields and the known target 
thicknesses, we determine the ratio of cross sections for diffractive dissociation 
on platinum and carbon (the two targets were subjected to essentially the 
same beam flux). The exponents a are then calculated using the cross section 
dependence a oc A a . The results are listed in Table 2 and compared with CT 
theoretical predictions [11]. The results are consistent with those expected 
from color-transparency calculations and clearly inconsistent with a values 
for incoherent scattering observed in other hadronic interactions. 



The Photon Light-Cone Wave Function 

The photon I.e. wave function can be described in a way similar to that of the 
pion except that it has two major components: the electromagnetic and the 
hadronic. Being a gauge field capable of piont-like coupling it has also a point 
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Fig. 5. Diagrams for photon dissociation 

bare-photon component. Consequently, the photon light-cone wave function 
in the pQCD regime can be expanded in terms of Fock states: 

■0 7 = a\jp) + b\l + l~) + c\l + l~^) + {other e.m.) 

+d\qq) + e\qqg) + (other hadronic) + .... (4) 

where \^ p ) describes the point bare-photon and |/ + /~) stands for |e + e~), |yU, + /z~) 
etc. The photon can be described as containing \qq) components even in the 
non-pQCD regime of the vector-meson dominance model [16], in particular the 
p, u,<f) and J/ip mesons. There is also a \qq) component with large relative 
transverse momentum that is connected to the point-like (direct) photon. The 
wave function of the photon is very rich: it can be studied for real photons, for 
virtual photons of various virtualities, for transverse and longitudinal photons 
and the hadronic component may be decomposed according to the quarks fla- 
vor. Strange quarks will reflect the strange content of the photon and charmed 
quarks will reflect the smaller sized charmed components. 

The point-bare photon does not have internal structure, it can only Compton- 
scatter and will not be studied here. The other electromagnetic Fock states 
begin with and continue to more complex systems. The interaction of 

these Fock states with the target is expected to be purely electromagnetic, as 
shown in figure 5(a). These were studied extensively through measurements of 
Bethe-Heitler leptoproduction [17] and of F^qed [18] • Exclusive measurements 
of the |/ + /~) wave function will complement these inclusive measurements, may 
provide more tests of QED and will help define the analysis tools for the more 
complex hadronic structure. 

The hadronic components interact according to the diagram of Fig. 5(b). The 
distribution amplitudes for these components, integrated over kt were calcu- 
lated by Petrov et al. [19] using the instanton model and are shown in figure 6. 
Balitsky et al. [20] predict for real photons a distribution amplitude identical 
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Fig. 6. The Photon wave functions, Petrov et al. Real photons (Solid line), virtual 
photons Q 2 = 250 MeV 2 (dashed line), Q 2 = 500 MeV 2 (dotted line). 

to the pion asymptotic function [2-4] (see Fig. 1). 

While a variety of cross sections, form factors etc. depend on the light-cone 
wave function they are usually not sensitive to its structure. This is because it 
normally enters the calculations in an integrated form. Such are also the pho- 
ton inclusive structure functions [18]. While this is an advantage when one 
needs calculations that do not depend on the internal structure of the photon, 
it is a disadvantage if we want to understand it. A differential measurement 
which will be sensitive to the u and k t dependence of <fi{u, Q 2 ) will test the 
most fundamental description of the photon internal structure. Such measure- 
ments will make it possible to compare with the theoretical predictions and 
determine the regime of their validity. 

The experimental program is presently in preparation and is based on diffrac- 
tive dissociation of a real or virtual photon into di-leptons or di-jets. The 
electromagnetic component will be studied through measurements of the u 
distribution for pure diffractive (e + e~) elastic photoproduction and the 

kt distribution of these events. Similar measurements of diffractive hadronic 
di-jets will be used to study the hadronic component. 

I would like to acknowledge the efforts of the E791 collaboration, of which 
I am a member, for the data presented in this work and the data analysis 
performed by my graduate student R. Weiss-Babai. 
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